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Two recent experiments have reported the trapping of individual atoms inside optical resonators 
by the mechanical forces associated with single photons [Hood et al.. Science 287, 1447 (2000), 
Pinkse et al., Nature 404, 365 (2000)]. Here we analyze the trapping dynamics in these settings, 
focusing on two points of interest. Firstly, we investigate the extent to which light-induced forces 
in these experiments are distinct from their free-space counterparts, and whether or not there are 
qualitatively different effects of optical forces at the single-photon level within the setting of cavity 
QED. Secondly, we explore the quantitative features of the resulting atomic motion, and how these 
dynamics are mapped onto experimentally observable variations of the intracavity field. Towards 
these ends, we present results from extensive numerical simulations of the relevant forces and their 
fluctuations, as well as a detailed derivation of our numerical simulation method, based on the full 
quantum-mechanical master equation. Not surprisingly, qualitatively distinct atomic dynamics arise 
as the coupling and dissipative rates are varied. For the experiment of Hood et al., we show that 
atomic motion is largely conservative and is predominantly in radial orbits transverse to the cavity 
axis. A comparison with the free-space theory demonstrates that the fluctuations of the dipole force 
are suppressed by an order of magnitude. This effect is based upon the Jaynes-Cummings eigenstates 
of the atom-cavity system and represents distinct physics for optical forces at the single-photon level 
within the context of cavity QED. By contrast, even in a regime of strong coupling in the experiment 
of Pinkse et al., there are only small quantitative distinctions between the potentials and heating 
rates in the free-space theory and the quantum theory, so it is not clear that a description of this 
experiment as a novel single-quantum trapping effect is necessary. The atomic motion is strongly 
diffusive, leading to an average localization time comparable to the time for an atom to transit freely 
through the cavity, and to a reduction in the ability to infer aspects of the atomic motion from the 
intracavity photon number. 
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I. INTRODUCTION 



An exciting advance in recent years has been the in- 
creasing ability to observe and manipulate the dynami- 
cal processes of individual quantum systems. In this en- 
deavor, an important physical system has been a single 
atom strongly coupled to the electromagnetic field of a 
high-Q (optical or microwave) cavity within the setting 
of cavity quantum electrodynamics (cavity QED). [y^] 
Here the coupling frequency of one atom to a single mode 
of an optical resonator is denoted by go (i-e., 2(7o is the 
one-photon Rabi frequency), with the regime of strong 
coupling defined by the requirement that ^ (7,^), 
where 7 is the atomic decay rate to modes other than the 
cavity mode and k is the decay rate of the cavity mode 
itself. In this circumstance, the number of photons re- 
quired to saturate an intracavity atom is no ~ 7^/5o 1 
and the number of atoms required to have an appreciable 
effect on the intracavity field is Nq ^ k-I / Qq 1 ||]- 

Although there have been numerous laboratory ad- 
vances which demonstrate the effect of strong coupling 
on the internal degrees of freedom of an atomic dipole 
coupled to the quantized cavity field (i.e., ^ k, 7), the 
consequences of strong coupling for the external, atomic 
center-of-mass motion with kinetic energy Ek have only 
recently been explored experimentally ft-S] . In a regime 



of strong coupling for the external degrees of freedom, 
go> Ek/h, a single quantum is sufficient to profoundly 
alter the atomic center-of-mass (CM) motion, as an atom 
moves through a region of spatially varying coupling co- 
efficient g{r) = goipir) [e.g., as arises in the Gaussian 
mode of a Fabry-Perot cavity, 'il'{r)]. 



Perhaps most strikingly, the spatial variation of the 
cavity mode can lead to a confining potential sufficient 
to trap an atom within the cavity mode even for a sin- 
gle quantum of excitation of the atom-cavity system, as 
first discussed in the work of Refs. [|io|,0. This is il- 
lustrated in Fig. 1, which shows the the possibility for 
trapping by excitation to the lower component |— ) in the 
Jaynes-Cummings manifold of eigenstates. Modifications 
of the atomic CM dynamics can in turn significantly al- 
ter the cavity field. This situation is very different from 
the usual case for trapped atoms or ions in fixed exter- 
nal potentials, in that here the confining field and the 
atomic motion can be strongly interacting, in which case 
the overall state of the system must be determined in a 
self-consistent fashion. 
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FIG. 1. The energy-level diagram for the coupled 
atom-cavity system, as a function of the atom's radial po- 
sition p. When the atom is near the cavity center, driving at 
frequency tjp populates the state |— ) to trap the atom. Here 



'^(p,c,a) ^(probe, cavity , atom) 



of the text. 



The experimental requirements to investigate strong 
coupling for both the internal and external degrees of 
freedom are stringent [namely, g > {Ek/h,j,K)], and 
have required the integration of the techniques of laser 
cooling and trapping with those of cavity QED, as was 
initially achieved in 1996 (l2j and as illustrated in Fig. 
2. Mechanical effects due to strong coupling with single 
quanta were first observed in 1998 [Q, in an experiment 
with peak coupling energy hgo ~ 5 mK and with initial 
atomic kinetic energy — 400 fi K. 

Following this theme, two groups recently reported 
trapping of single atoms with intracavity fields at the 
single-photon level, beginning with the work of Ref. H] 
and culminating in that of Refs. That such trapping 

might be possible in these experiments is indicated by the 
fact that the ratio R of initial atomic kinetic energy Ek 
to the coherent coupling energy hgo^R = Ek/hgo, is less 
than unity. For the work in Refs. R — 0.06, while 

for that in Ref. i? ~ 0.27. Although these ratios are 
indicative of the possibility of trapping with single quanta 
in cavity QED, the actual forces and confining potentials 
are somewhat more complex to analyze, as we shall see. 
Moreover, beyond providing single-quantum forces suffi- 
cient for atomic localization, strong coupling also means 
that the presence of one atom can significantly modify 
the intracavity field, thereby providing a means to track 
atomic motion by way of the light emerging from the 
cavity. 

To understand the basic scheme for trapping of single 
atoms with single quanta in cavity QED, consider the 
energies HP± for the first excited states |±) of the atom- 
cavity system. Along the radial direction p — ^ ip- + 
and for optimal x (standing- wave) position, P±{p) has 



the spatial dependence indicated in Fig. 1, which ne- 
glects dissipation. The ground state of the atom-cavity 
system is |a,0); the atom is in its ground state a, and 
there are no photons in the cavity. For weak coupling 
(atom far from the cavity mode center), the first two 
excited states are that of one photon in the cavity and 
the atom in the ground state, \a, 1), and of the atom in 
the excited state e with no photons in the cavity, |e,0). 
These two states are separated by an energy hAac, where 
Aac = cavity— ^ atom IS the detuning between the "bare" 
(uncoupled) atom and cavity resonances. 

FIG. 2. Experimental schematic for the case of Hood et al. 
Atoms are captured in a magneto-optical trap (MOT), and 
dropped or launched through a high-finesse optical cavity. A 
single atom (trace with arrow) transiting the cavity mode al- 
ters the measured transmission of a probe beam through the 
cavity. In the experiment of Pinkse et al, rubidium atoms 
are captured in a MOT below the cavity and launched up- 
ward through it. 

As an atom enters the cavity along p it encounters 
the spatially varying mode of the cavity field, and hence 
a spatially varying interaction energy hg{r), given by 
g{r) = 5oCOs(A;a;/)exp(-(?/2 + z^)/wl) {k = 2it/\). The 
bare states map via this coupling to the dressed states 
|±) shown in the figure, with energies 



/3± 



^atom ~(~ ^cavity 



ffacll/2 



(1) 



Our interest is in the state |— ); the spatial dependence of 
the energy hjd- (r) represents a pseudopotential well that 
can be selectively populated by our choice of driving field 
Sprobeit) and Aprobe to trap the atom, as first suggested 
by Parkins p^ . The system is monitored with a weak 
probe beam as an atom enters the cavity mode; detection 
of an atom transit signal triggers an increase in driving 
strength to populate the state |— ) and trap the atom. Be- 
cause the experiments in the optical domain have atomic 
and cavity decay times (k~^,7~^) that are small com- 
pared to the time r for motion through the cavity field, 
the atom-cavity system must be continually re-excited by 
way of Sprobe, thereby providing an effective pseudopo- 
tential on time scales St such that {k^^, 7^^)^ St ^ r. 

Although a full theory based on the preceding discus- 
sion is sujficient to provide detailed agreement with the 
experimental observations of Refs. (^J^,^ (as we shall 
show in subsequent sections), it is reasonable to ask to 
what extent such a theory based on the interactions in 
cavity QED is necessary. In particular, it might well 
be that the well-established theory of laser cooling and 
trapping in free space jlj] could provide an adequate de- 
scription of the potentials and heating rates, with the 
cavity merely providing a convenient means for attaining 
a strong drive field. With respect to the experimental re- 
sults of Pinkse et al. (Ref. |^), we find that this is in fact 
largely the case; there are only small quantitative distinc- 
tions between the free-space theory and the appropriate 
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quantum theory. One interesting feature to note in this 
experiment is enhanced coohng of the atomic motion rel- 
ative to the parameters of Hood et al. This effect, 
which enables trapping in this parameter regime, arises 
through cavity-mediated cooling For these pa- 

rameters, the average localization time from simulations 
is extended by 75% relative to the equivalent free-atom 
signal; both these times are shorter than the time for an 
atom to transit freely through the cavity. 

By contrast, in the regime of the experiment of Hood 
et al. (Ref. the cavity QED interactions result in a 
strong suppression of dipole heating along the cavity axis 
relative to the free-space theory, which has a strong ef- 
fect on both the duration and character of the observed 
atom transits. In the cavity QED setting it becomes 
possible to create a potential deep enough to trap an 
atom without simultaneously introducing heating rates 
that cause rapid escape from that potential. For these 
parameters, the average experimentally observed local- 
ization time is a factor of 3.5 longer than the equivalent 
free-atom average. The results of extensive numerical 
simulations of trapping times and radial oscillation fre- 
quencies, and their validation by way of comparisons to 
experimentally measured distributions, demonstrate the 
essential role of the single-photon trapping mechanism in 
the experiment of Ref. Q . At root is the distinction be- 
tween the nonlinear response of an atom in free space and 
one strongly coupled to an optical cavity. For these ex- 
perimental parameters, the eigenvalue structure of Fig. 1 
leads to profound differences between the standard the- 
ory of laser cooling and trapping, and the extension of 
this theory to the regime of strong coupling in cavity 
QED. 

Note that prior experiments in our group have con- 
firmed that the full quantum treatment of the one-atom 
master equation in cavity QED is required for a descrip- 
tion of the dynamics associated with the internal degrees 
of freedom for a single atom in an optical cavity in the 
regime g > (7, k). These experimental confirmations 
come by way of measurements of the nonlinear suscep- 
tibility for the coupled system in settings close to that 
for the experiment of Ref. [pp|Jl7[|. A principal goal 
of this paper is to investigate the extent to which a the- 
ory of atomic motion within the setting of cavity QED is 
likewise a necessary component in describing the center- 
of-mass dynamics for the experiments of Ref s. 

A second goal is to examine the related question of the 
extent to which inferences about atomic motion within 
the cavity can be drawn from real-time observations of 
the cavity field, either via photon counting or hetero- 
dyne detection of the cavity output. The interactions 
in cavity QED bring an in principle enhancement in the 
ability to sense atomic motion beyond that which is oth- 
erwise possible in free space. Stated more quantitatively, 
the ability to sense atomic motion within an optical cav- 
ity by way of the transmitted field can be character- 
ized by the optical information 

which, roughly speaking, is the maximum possible num- 



ber of photons that can be collected as signal in time At 
with efficiency a as an atom transits between a region 
of optimal coupling go and one with g{r) <Si On. A key 
enabling aspect of the experiments in Refs. is that 
TZ —g^/n ^ ('t,7), leading to information about atomic 
motion at a rate that far exceeds that from either cav- 
ity or spontaneous decay (as in fluorescence imaging). 
In practice, for detection strategies employed experimen- 
tally, information is extracted at a somewhat lower rate. 
For example, in the experiment of Hood et al. 0, the 
photon count rate would be (2.7 x 10^/s) (including the 
overall escape and detection efficiency a « 0.15), while 
for the experiment of Pinkse et al. |^ it is (2.2 x 10^/s) 
(including an estimated overall escape and detection effi- 
ciency a « 0.11) For time scales At ^ lOfi s as rele- 
vant to the following discussion, atomic motion through 
the spatially varying cavity mode leads to variations in 
the transmitted field that can be recorded with a high 
signal-to-noise ratio, namely, a signal of 2.7 x 10^ pho- 
tons for the experiment of Hood et al. and 2.2 x 10^ 
for that of Pinkse et al., where each is calculated for an 
intracavity field strength of one photon. 

The value of the optical information itself does not tell 
the complete story. For cavity QED experiments like 
those considered here, one records either the sequence 
of photoelectric counts or the heterodyne current ver- 
sus time, from which necessarily only limited inferences 
about atomic motion can be drawn. However, if center- 
of-mass dynamics (i.e., axial and radial motions) occur 
on well-separated time scales, then it is reasonable to sug- 
gest that appropriate signal processing techniques could 
extract information about these motions from the single 
time sequence of the photocurrent i{t). Such process- 
ing could presumably occur in real time if aTZ is much 
faster than the rates for radial and axial motion [e.g., the 
oscillation frequencies {,fr,fa) in a potential well, with 
fr <SC fa\- Unfortunately, in neither experiment is 
alZ large enough to resolve the axial dynamics directly, 
so the task of disentangling the radial and axial motion 
signals becomes more difficult, and theoretical simula- 
tions of the experiment become useful in understanding 
the nature of the observed transmission signals. 

This difficulty arises in the experimental regime of 
Pinkse et al. For these parameters, axial heating 
leads to frequent bursts of large-amplitude motion along 
the cavity axis, with envelopes extending over time scales 
comparable to those for radial motion. Consequently, at 
experimental bandwidths (averaging times), both types 
of motion give rise to qualitatively similar modulations 
in the measured transmission signal. Furthermore, mo- 
tion in the radial direction has a strong diffusive compo- 
nent, giving rise to a wide spread of time scales for ra- 
dial motion. Our simulations discussed in Sec. V suggest 
that for these parameters, short-time-scale modulations 
(< 300 /zs) tend to be mostly due to bandwidth averag- 
ing over axial motion, while longer (> 500 /i s) variations 
such as presented in Fig. 2 of Ref. |^ typically reflect ra- 
dial motion, though these long-time-scale variations are 
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generally modified in amplitude by the presence of axial 
motion. Modulations on intermediate time scales appear 
ambiguous in their dynamical origin. 

By contrast, as shown in Ref. |^], for the parameters of 
Hood et al. atoms are well localized along the standing- 
wave direction throughout most of the trapping interval, 
with axial motion giving rise to negligible signal until 
finally rapid axial heating leads to atomic escape. Con- 
sequently, observed variations in the photocurrent i{t) 
are simpler than those of Ref. |^ , and directly yield the 
radial atomic position. Furthermore, in this experiment 
the radial oscillation frequency is large compared to the 
spontaneous emission heating rate, meaning that the re- 
sulting atomic motion is largely conservative (rather than 
diffusive) in nature, taking place in a known potential (as 
demonstrated both experimentally and by way of numer- 
ical simulation). Hence, from i{t) it becomes possible to 
make detailed inferences about the radial motion, even 
to the point of real-time observations of the anharmonic 
motion of a single atom and of the reconstruction of ac- 
tual atomic trajectories. 

The structure of the paper is as follows. Following this 
introduction, in Sec. H we present a detailed description 
of our theoretical model and its use for the implementa- 
tion of numerical simulations. Section HI compares ef- 
fective potentials and momentum diffusion rates derived 
for the two experiments, along with their analogs for the 
hypothetical case of an equal-intensity free-space trap. 
These calculations explore the distinction between quan- 
tum and classical, and also give insight into the nature 
of atomic motion expected in both experiments. Sample 
simulated trajectories are presented for both cases. In 
Sec. IV we present experimental and simulation results 
for the case of Hood et al., which serve both to verify the 
simulations and also to demonstrate important features 
of the resulting motion. Sec. V gives the application of 
the same tools to analyze the experiment of Pinkse et 
al.; we see that standing- wave motion and diffusive ra- 
dial motion complicate the correlation between atomic 
position and detected field in this case. Finally, axial 
motion is explored in more depth, and Fourier analysis 
of our simulations show that oscillations of comparable 
amplitude and frequency should be visible for both atoms 
confined (but heated) within a well, and atoms skipping 
along the standing wave. 



Principal findings 

The theoretical treatment and numerical simulation of 
the motion of a single atom strongly coupled to an op- 
tical cavity, as described in Sec. II, lead to a surpris- 
ingly rich range of often qualitatively different dynamics. 
The motion may be essentially conservative and tightly 
confined around antinodes of the standing wave, or es- 
sentially dissipative and diffusive and involve interesting 
flights between different potential wells of the standing 



wave. Indeed we find that the existing experimental re- 
sults of Hood et al. and Pinkse et al. exemplify these very 
different dynamical regimes. Key features of the atomic 
motion in both experimental regimes are addressed as 
follows: 

Figures 3-5 and their associated discussion in Sec. Ill 
elucidate the nature of the trapping potential and mo- 
mentum diffusion in an optical cavity as opposed to a 
free space standing wave. In particular we find that, 
even when the atom-cavity system is strongly coupled 
and driven such that it has a mean intracavity photon 
number of roughly 1, the trapping potential and momen- 
tum diffusion may be only slightly different from those in 
a free-space standing wave, and in fact this is the case for 
the parameters of Pinkse et al. On the other hand, for 
the parameters of Hood et al. the usual fluctuations of 
the dipole force along the standing wave are suppressed 
by an order of magnitude, which to our knowledge rep- 
resents qualitatively new physics for optical forces at the 
single-photon level within the context of cavity QED. We 
show that in the parameter regime of Pinkse et al. the 
heating rates are such that the atom could be expected 
to gain energy equal to a signiflcant fraction of the to- 
tal trapping potential during a single motional oscillation 
period for both axial and radial motion. By this measure 
the heating rates in the experiment of Hood et al. are 
much slower, indicating more nearly conservative motion, 
and this could be expected to have a profound effect on 
the qualitative nature of the dynamics in the two exper- 
iments. 

Figures 6 and 7 and the corresponding text in Sec. Ill 
present simulated transits for both experiments, and dis- 
cuss the qualitative features of atomic dynamics in both 
cases. For the parameter regime of Hood al., conserva- 
tive radial motion dominates diffusion and standing-wave 
motion, with atomic trajectories localized at peaks of a 
single standing-wave antinode. Atoms trapped with the 
mean trapping time execute several radial orbits. The 
eventual escape is typically due to heating along the cav- 
ity axis. By contrast, for the experiment of Pinkse et al., 
a trajectory of typical duration, as in Fig. 7(a), does not 
experience a complete radial orbit and in fact resembles 
a scattering event, with a large contribution from radial 
diffusion as well. For these events the observed localiza- 
tion time is comparable to the time for free flight through 
the cavity. Axially the simulations show that in longer 
duration transits the atom frequently skips between wells 
of the standing-wave potential due to repeated heating 
and recooling. 

Section IV, with Figs. 8-10, presents a more detailed 
and quantitative investigation of trapping and motional 
dynamics for the experiment of Hood et al. The ability 
of our simulations to closely reproduce the mean trap- 
ping times observed in the experiment provides evidence 
of their accuracy and utility. As illustrated in Fig. 9, 
the triggering strategy leads to signiflcant modiflcations 
of the distribution of residence times within the cavity. 
The essentially conservative nature of the dynamics and 
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the strong axial confinement make it possible to confi- 
dently ascribe oscillations in the transmitted intensity to 
radial motion of the atom. As shown in Fig. 10, the 
experimentally observed oscillations are consistent with 
the calculated potential. The conservative nature of the 
motion is further confirmed by the separation of orbital 
periods by angular momentum that is also apparent in 
this figure. 

Section V, with Figs. 11-15, presents a detailed analy- 
sis of trapping and motional dynamics for the experiment 
of Pinkse et al. Again, our simulations are sufficient 
to reproduce the reported mean localization time. In 
this case, the triggering strategy leads to relatively mi- 
nor modifications of the distribution of residence times 
for an atom within the cavity. In this case the dissipa- 
tive nature of the evolution is significant; essentially no 
long-term localization is observed if the sign of the fric- 
tion coefficient is reversed, whereas this has little effect 
in the parameter regime of Hood et al. These largely 
dissipative and diffusive motional dynamics are found to 
have significant effect on the information about the mo- 
tion that is available in the transmitted field. For those 
events with a long localization time, the axial motion of 
the atom is repeatedly heated and cooled, resulting in 
slow variations in envelope of the amplitude of the rapid 
oscillations of the transmitted light. The time scale of 
these variations is comparable to that for radial motion 
of the atom. There are thus no unambiguous signatures 
for radial motion and only longer time scale excursions of 
the atom in the radial potential lead to variations of the 
output field that may be confidently ascribed to the ra- 
dial motion. Likewise, although information about axial 
motion is also available in the output light, we find that 
it is in general difficult to distinguish large oscillations in 
a single well of the axial potential from free flight over 
several wells as attempted in Ref . [|j . 



II. THEORETICAL MODEL AND NUMERICAL 
SIMULATIONS 

In this section we outline the derivation from the full 
quantum-mechanical master equation of the "semiclas- 
sical model" for the atomic motion used in Ref. 0]. It 
turns out that this model is able to reproduce the ex- 
perimental observations very accurately. Note that here 
the term "semiclassical" refers to approximations with 
respect to the atomic center-of-mass motion, and not to 
the internal degrees of freedom, for which the full quan- 
tum character is retained. This situation should not be 
confused with the semiclassical theory of cavity CJED 
for which expectation values of field operators O field 
and atomic operators O field are assumed to factorize, 

{OfieidOatom) = {O field) {O atom) ] no such approxima- 
tion is made here. To distinguish these two cases, we 
introduce the term quasiclassical for the case of atomic 
motion. 



The validity of the quasiclassical model depends on a 
separation of time scales between the atomic motion and 
the cavity and internal atomic dynamics. We adapt the 
work of Dalibard and Cohen- Tannoudji jlj] to the situ- 
ation of a quantized cavity mode. A similar derivation 
in the bad-cavity limit appears in |ig| ]. The details of 
the derivation are essentially unchanged from free space, 
since the terms of the master equation which refer to the 
dynamics of the cavity have no explicit dependence on 
the operators describing the atomic motion. However, 
we do find conditions for the validity of the approxima- 
tion for this system which depend on the properties of 
the cavity. Finally, we describe in more detail the nu- 
merical simulations of the resulting model first presented 
in Ref. |71|. These simulations are of the kind discussed 
in Refs. ||2^,|ltl. 

An analytical calculation of force, momentum diffu- 
sion and friction coefficients for the quasiclassical model 
of atomic motion in the low driving limit was derived 
by Horak and co-workers |l^j2^, who found a regime in 
which the steady-state temperature scaled as the cav- 
ity decay rate. This allows cooling of the atom below 
the Doppler limit, so long as the cavity can be made to 
have lower loss than the atom. However, the parame- 
ters of Refs. 1^,^ are very far from this low driving limit. 
Hence we employ numerical techniques based on solv- 
ing the appropriate master equations by expansions in 
terms of Fock states of the cavity field P3] . Note that a 
very early contribution developed a different theoretical 
framework and numerical scheme for calculating the force 
and friction (but not the momentum diffusion) of an atom 
in a cavity (or "colored vacuum" ) |^,^ . Very recently, 
Vuletic and Chu [l6| found cavity-mediated cooling in a 
slightly different regime to that considered by Horak et 
al. 



A. Model of atom-light interaction in a cavity 

The Hamiltonian for a two-level atom interacting with 
a single mode of the electromagnetic field in an optical 
cavity using the electric dipole and rotating-wave approx- 
imations (in the interaction picture with respect to the 
laser frequency) is 

+hg{r){a^(j + (j''a) + fi{£a^ +£*a) . (2) 

This is the familiar Jaynes-Cummings Hamiltonian mod- 
ified to take into account the external degrees of freedom 
of the atom and the spatial variation of the cavity mode. 
The first term is the kinetic energy of the atom, and 
the next two terms are the energy in the internal state 
of the atom and the cavity excitation. The fourth term 
describes the position-dependent interaction of the cav- 
ity mode and the atomic dipole. It is important to note 
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that r and p are operators. Thus, for example, the ex- 
act strength of the couphng between the atomic internal 
state and the cavity field depends on the shape of the 
atomic wave packet, which is in turn determined by the 
mechanical effects of the cavity field. Some implications 
of this Hamiltonian are considered in detail by Vernooy 
and Kimble p6t . The Hamiltonian has been written in 
terms of cavity and dipole operators that rotate at the 
frequency of the probe field ujprote ■ The real atomic tran- 
sition (cesium in Ref . and rubidium in Ref . |§1 ) in fact 
involves several degenerate magnetic sublevels, but we 
assume that the cavity is driven by circularly polarized 
light and that the atom is optically pumped such that it 
occupies an effective two-level system described by the 
dipole operator cr with the quantization axis along x. 

Dissipation in the system is due to cavity losses and 
spontaneous emission. By treating modes external to 
the cavity as a heat reservoir at zero temperature in the 
Born, Markov, and rotating-wave approximations, it is 
possible to derive the standard master equation for the 
density operator p of the system as 

= — [H, p] + K{2apa} — a) ap — pa) a) 



47r 



(3) 



The third and fourth terms describe the effect of spon- 
taneous emission on the atomic motion including the 
momentum kick experienced by the atom as a result 
of the spontaneous emission. The unit vector k is 
the direction of an emitted photon. The pattern of 
dipole radiation is accounted for by the angular factor 

5(k-x) = [l + (k-x)'l /2 pi. 



B. Quasiclassical motion of the center of mass 

It is possible to eliminate the internal and cavity dy- 
namics adiabatically in favor of the slower dynamics of 
the motional state in parameter regimes of direct rele- 
vance to current experiments. Intuitively, for the qua- 
siclassical approximation to work, the state of the atom 
needs to be sufficiently localized in position and momen- 
tum on the scales important to the problem so that it 
can be thought of as a classical particle. The conditions 
for adiabatically eliminating the internal and cavity dy- 
namics roughly correspond to this idea. It turns out that 
it is necessary first that exchanges of momentum with ei- 
ther the cavity field or by spontaneous emission into free 
space should result in momentum kicks that are small 
compared with the momentum spread Ap of the atomic 
Wigner function, thus 



For an atom which is in a minimum uncertainty state 
with respect to the position-momentum Heisenberg in- 
equality this requires that the state is localized to better 
than a wavelength. The atomic motional state will in 
general be a mixture allowing the position spread to be 
broader. However, this requirement means that the mo- 
tional state can be thought of as a probabilistic mixture 
of pure states localized to within a wavelength, and so 
places a limit on the coherence length of the motional 
state ||2^]. Second it is important that the range of 
Doppler shifts of the atom due to its momentum spread 
is small compared to the atomic and cavity linewidths, 
thus 



£2 ~ fcAp/m7 ~ fcAp/mK ^ 1. 



(5) 



In this paper it will be assumed that the root-mean- 
square atomic momentum obeys this inequality, thus 
making a low velocity approximation, but the arguments 
here can in fact be generalized to arbitrary mean veloc- 
ities of the atom [Q. The Heisenberg inequality means 
that this also requires a minimum position spread of the 
atom 



Ar 3> hk/TW-f, hk/mK. 



(6) 



These criteria are a simple generalization of the situation 
for laser cooling in free space which can be imagined as 
the situation k oo. The consistency of these condi- 
tions, which effectively put lower and upper limits on the 
atomic momentum spread, requires that 



h'-f 



« 1, 



< 1. 



(7) 



£i ~ h/Ap < 1. 



(4) 



The first of these conditions is well known for laser cool- 
ing in free space — the requirement that the recoil energy 
of the atomic transition be much lower than the Doppler 
energy, which effectively controls the limiting tempera- 
ture of the laser cooling. This condition is well satis- 
fied for heavy atoms such as cesium and rubidium and 
the optical transitions employed in cavity QED experi- 
ments considered here. The analogous condition brought 
about by the cavity dynamics requires that the recoil 
energy associated with exchanging excitation with the 
cavity field is much smaller than the energy width of 
the cavity resonance. Just as the first criterion implies 
that the atom still be in resonance with a driving field 
at its transition frequency after spontaneously emitting, 
the second criterion implies that absorbing or emitting 
a photon from the cavity will leave the atom near the 
cavity resonance. In the experiments of Refs. |^,0,D, 
K ~ 7, so that this second criterion does not place a 
stronger restriction on the validity of the approximations 
than the free-space limit. However, it is important to 
note that the design of the cavity, as well as the atom 
and transition that are chosen, now has an effect on 
the validity of the approximation. It would be possi- 
ble, for example, to change the cavity length in such a 
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way that the system moves from a regime in which the 
quasiclassical treatment is appropriate into one in which 
it is not. In practice for cold atoms cooled to roughly the 
Doppler limit (Ap2/2m ~ ^17, Hk) it will be the case that 
£1 ~ £2 ~ y/{h?k^ /2m) //17, -y/ {h?k^/2m) /hn, and so a 
consistent expansion should be to equal order in these 
small parameters. 

The derivation of Ref. may be applied to our 
problem, and proceeds by transforming the master equa- 
tion [Eq. ^ into an evolution equation for a Wigner op- 
erator, 



W{r,p,t) 



1 



1 



1 



d u{f-\ — mIpIt u) exp(— • u/h) 



(8) 



dt m df 



describing the complete state of the system. An approx- 
imate Fokker-Planck equation for the Wigner function 
describing the motional degrees of freedom alone is found 
by writing this equation as a Taylor expansion in terms 
of the small parameters ei and £2, and truncating that 
expansion at third order. The force operator is defined 
as the gradient of the atom-cavity coupling 



F{r) = -hgoVip{r){a)a + ah 



(9) 



It is possible to show that the Fokker-Planck equation 
for the atomic Wigner function / takes the form 



dpidp. 



dptdp- 



dpidr. 



(10) 

The quantities appearing in the Fokker-Planck equation can be calculated from the master equation for the internal 
and cavity degrees of freedom alone that is obtained by setting r to some real number value ro, and disregarding the 
kinetic-energy term. We define Ps{r) as the steady state of this master equation, with the steady-state expectation 
value of the arbitrary operator c given by (c)^^ =Tr(c/Os(r)). The parameters appearing in the Fokker-Planck equation 
can then be expressed as follows: 



(f) = Tr F(r'V,(r) 
'1 



dT 



{F,{T)F,iO)+F,{0)F,{T))^-^,<j,, 



r]ij = — / drr 



- (F, (r) F, (0) + F, (0) F, (r))^^ - 0,0, 



r., = ;^/ drr([F,(r),F,(0)]) 



Simple integrations give Exx — 2/5 and Eyy = 3/10 = 
Ezz and all other components of E are zero. Excepting 
the different definition of the force operator F, these are 
the expressions that can be derived in case of a free-space 
light field p^ . However, it is important to bear in mind 
the extra conditions on the validity of the adiabatic elimi- 
nation. The master equation [Eq. ^ means that the force 
expectation values and correlation functions can be very 
different from those that are calculated in free space. In 
practice, the contribution from the parametric tensor 77 is 
often smaller than that from the diffusion tensor D by a 
factor of order £, and is usually disregarded in treatments 
of free-space laser cooling Q . 

Thus, as assumed in earlier work, calculating the qua- 
siclassical motion of the atom in a cavity field only re- 
quires that the force and its correlation function be eval- 
uated for the full atom-cavity master equation. Such 
prior treatments assumed that the atom is motionless; 



however, they can be extended to atoms moving at some 
velocity under the same conditions • The diffusion 

coefficients may be found by first calculating the corre- 
lation functions via the quantum regression theorem and 
numerical integration, or directly via matrix-continued 
fraction techniques ]3^j3l|] . A matrix-continued frac- 
tion calculation requires that the field mode be periodic, 
and as such it only works along the standing-wave axis 
of the cavity mode. In directions perpendicular to this, 
the calculation of correlations from the master equation 
is essentially the only option if the atom is not slowly 
moving. 



C. Stochastic simulations of the quasiclassical model 

It is possible to recast the Fokker-Planck equation of 
Eq. nG into a simple set of stochastic equations which 
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describe atomic trajectories in the cavity field. These diffusion and friction tensors can be rewritten using the 
equations can be used to gain intuition about the atomic definition of the force operator [Eq. (^] 
motion and how it is affected by mechanical forces. The 



D = 



fc2 2 



dT 



- ($ [r] $ [0] + $ [0] $ [r])^^ - ($) 



ml J L 



m L J L J 



dTT{[^ [r],<i>[0]])^^ 



(11a) 



(lib) 



where $ = aV + a^a. Writing the parameters of the 
quasiclassical model in this form relies on the approxima- 
tion that the atom is slowly moving, namely, that it does 
not move a significant fraction of a wavelength during a 
cavity or atomic lifetime. Note that the functions ^ and 
X depend on position only through the coupling g = goip- 
They can be calculated efficiently by finding D^^ and T^x 
using matrix-continued fractions, and then dividing off 
the gradient factors. A matrix-continued fraction tech- 
nique cannot be used to find the other components of the 
momentum diffusion or the friction tensors directly, since 
the field mode is not periodic across the Gaussian profile 
of the mode. 

It is now straightforward to convert the Fokker-Planck 
equation for the Wigner function into an equivalent set of 
Ito stochastic differential equations. The resulting (Ito) 
equations are |Q 

dx — —pdt, (12a) 
m 

dp = -hgo{<i>)Vijdt - ■^x(f) (p- Vip) Vip 

+2hgoy/SXj^Vi;dWi + 2hk-/^J {a^a)VEdW. (12b) 

The Wiener increment dWi has the usual properties, in 
particular dW^ = dt. The vector dW is a vector of 
three such increments. The terms in the equation for the 
momentum are the mean radiative force, its first-order 
dependence on momentum, and its fluctuations due to 
the atom-cavity system and due to the coupling to free 
space, respectively. These equations depend on the quan- 
tities ($),XjC and (crV), which are functions of position 
through g only. A straightforward simulation of these 
equations only needs to store ordered look-up tables of 
these quantities for given values of g, rather than for all 
possible values of r. All of the other quantities that ap- 
pear, including g, are simple functions of r and p. At each 
time step the algorithm searches the look-up table for the 
current value of starting from the previous value, and 
reads off the current values of (?'),X,'Cj (f^o")- A linear 
interpolation for the two closest values of g was used, but 
more sophisticated interpolation schemes could be imple- 
mented. Since g will not change by a large amount in any 



one timestep the search can be very efficient; a routine 
from Ref. |33j was used for this. In the low- velocity limit 
of the quasiclassical theory, these stochastic differential 
equations describe all the motional dynamics of the atom 
inside the cavity. The term proportional to rjij leads to 
correlations between the atomic position and momentum. 
The effect of rjij is typically small compared to friction 
and diffusion and has been ignored for the moment as 
is common practice in free-space standing waves. Terms 
in the SDE corresponding to the 77 term in the Fokker- 
Planck equation could easily be added. This would mean 
adding a new noise source which would affect the evolu- 
tion of the position as well as the momentum. 



III. APPLICATION OF THE MODEL TO 
EXPERIMENTAL REGIMES 

A. Potentials and heating rates for atomic motion 

The "quasiclassical" model discussed in the previous 
Sec. II can give us a great deal of information about the 
nature of the dynamics that may be expected in the pa- 
rameter regimes relevant to the experiments of Hood et 
al. and Pinkse et al. |^ . In particular we are interested 
in whether quantization of the cavity field leads to any 
significant change in the dynamics, in the sense of asking 
whether the atomic motion is very different in the cav- 
ity from what it would be in a free-space standing wave 
of the same intensity and geometry as the cavity mode. 
Second, we can investigate the nature of the resulting 
atomic motion in the cavity field, which can be either 
predominantly conservative or significantly diffusive and 
dissipative, depending on the particular parameters of 
interest. 

To obtain a feel for the type of atom dynamics ex- 
pected, effective potentials and heating rates were calcu- 
lated for both axial and radial directions of motion. The 
effective potential of the atom in the cavity field may be 
calculated from the force by 

U{r) = - [ F{r) ■ dr. 
Jo 
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The heating rates represent the average increase in the 
motional energy due to the momentum diffusion at a 
given position f, and may be calculated from the dif- 
fusion tensor according to 



dE 
'dt 



(r) = Tr[D(f)]/r 



Thus the axial potential at the center of the mode is 
U (0, x) = — F{0, x')dx' and the associated axial heat- 
ing rate is dE{0, x)/dt — DxxiO,x)/m. These quan- 
tities along with their radial equivalents U{p,0) and 
dE{p,0) /dt are plotted in Fig. 3 for the parameters of 
Hood et al. [0. The force and momentum diffusion co- 
efficient for the cavity system were calculated according 
to the formulas described above by numerical techniques 
based on Ref. The field state is expanded in terms 
of number states, and truncated at an appropriate level 
and a matrix-continued fraction algorithm is used to cal- 
culate D. The axial potentials and heating rates have 
A/2 = 426 nm periodicity inherited from the standing- 
wave field strength. Observe that the axial heating rates 
have minima at both field antinodes and field nodes. 
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FIG. 3. Effective potentials t/e// and heating rates dE/dt 
in the radial and axial directions for the experiment of 
Hood et al. (solid traces). The cavity field has a Gaus- 
sian waist W(i — 14 in the radial direction. The ax- 
ial standing wave has antinodes at x — (0, ±0.426)/im and 
nodes at s = ±0.213/im. All quantities are calculated for 
Ap,.obe/27r = -125 MHz and Aac/^-K = -47 MHZ, with an 
empty cavity mean-field strength of m = 0.3 photons. For 
comparison, corresponding quantities for an equivalent clas- 
sical free-space trap are shown as dashed traces. Note that 
the axial heating in the cavity trap is tenfold smaller, greatly 
enhancing the trap lifetime. 

The first thing to note is that the axial and radial 
heating rates are very different. In the radial direction, 
heating is dominated by diffusion due to spontaneous- 
emission recoils. Axially, however, the reactive or dipole 
fluctuation component of the diffusion dominates. This 



is because the reactive component is proportional to the 
gradient of the field squared, which is much larger for the 
axial direction where variations are greater (by a factor 
of 27niio/A). This contribution also has the property that 
it does not saturate with the atomic response. 

It is already clear that it should be possible to trap in- 
dividual atoms, since the potential depth of roughly 2.5 
mK is greater than the initial energy of the atoms in the 
experiment (around 0.46 mK) and the heating rate in 
the radial potential is relatively slow. Over 50 fis (a time 
scale over which the atomic motion is strongly affected 
by the potential) the total heating will typically still be 
small compared to the depth of the potential. However, 
the importance of the quantum character of the relevant 
fields or phenomena is not ensured by the statement that 
trapping occurs with a mean field strength of m ~ 1 
photon, since this is trivially the case in an equivalent 
free-space volume for a field of the same intensity as that 
inside the cavity. 
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FIG. 4. Effective potentials Ueff and heating rates dE/dt 
in the radial and axial directions for the experiment of Pinkse 
et al. (solid traces). The cavity field has a Gaussian waist 
Wo = 29 pm in the radial direction. The axial standing wave 
has antinodes at a; = (0, ±0.390) pm and nodes ai x — ±0.195 
fim. All quantities are calculated for Aprobe/'2'K = —45 MHz 
and Aac/2TT — —40 MHz, with an empty cavity photon num- 
ber n = 0.9. For comparison, corresponding quantities for 
an equivalent classical free-space trap are shown as dashed 
traces. Note that the potential depths and heating rates are 
comparable in the cavity QED and free-space cases. 

In order to see whether a full quantum description of 
the atom-cavity is necessary in order explain observed 
effects. Fig. 3 also shows the values calculated for an 
atom in an equivalent free-space standing wave, calcu- 
lated by standard techniques |2|] . This free-space stand- 
ing wave has the same geometry as the cavity mode, and 
the same peak field strength go|(a)P(0, 0). The detun- 
ing between the free-space field and the atom is chosen 
to be Aprobe- Pcrhaps surprisingly, the only large differ- 
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ence between the two models is in the axial heating rate, 
where a strong suppression of the axial heating is seen 
in the quantum calculation. This suppression is an effect 
of the quantized nature of the intracavity field. The self- 
consistent coupling of the cavity field and atomic position 
(in a semiclassical sense) cannot explain this suppression; 
in fact, by itself this coupling would lead to an increase in 
diffusion over the free-space case, since the atomic mo- 
tion within the cavity induces steeper gradients in the 
field. The suppression of diffusion is then evidence that 
it is necessary to use a fully quantum description, and 
speak of single photons rather than classical fields for 
these experimental parameters. As discussed in Ref. [Q, 
this suppression of the axial heating was essential for the 
trapping of atoms in the cavity. Thus for these exper- 
imental parameters, the eigenvalue structure of Fig. 1 
leads to profound differences between the standard the- 
ory of laser cooling and trapping and the extension of this 
theory to the regime of strong coupling in cavity QED. 
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FIG. 5. Comparison of effective potentials and heating 
rates in the experiments of Hood et al. (solid line) and Pinkse 
et al. (dash-dotted line). Heating rates are shown in units of 
trap depths per harmonic oscillation period (in the appropri- 
ate trap dimension), providing a direct measure of the degree 
to which oscillatory motion can be expected to be conserva- 
tive in nature. Note that differences in wo and A between 
the two experiments lead to quite different radial widths and 
slightly different axial periodicities for the quantities plotted. 



By way of comparison, the same quantities are plotted 
for the parameters relevant to Pinkse et al. ^ in Fig. 4. 
[Q. The smaller value of go in this experiment leads to 
a smaller effective potential, since the spatial gradients 
of the dressed state energy levels (which lead to the po- 
tential) are proportional to go. More importantly, the 
diffusion values calculated from the full quantum model 
discussed above are now little different from those of the 
equivalent free-space standing wave. This lack of a clear 
difference in potentials or diffusion indicates that the 
quantized nature of the field is not required to explain 



the radial trapping observed in Ref. Note that the 
resulting axial heating rates are essentially the same as 
those of Ref. in absolute magnitude; however, in Ref. 
[Q the potential was made deeper without the expected 
corresponding increase in diffusion. For the parameters 
of Ref. 1^ one additional interesting feature appears — 
enhanced cooling of the atom motion relative to the pa- 
rameters of Ref. . This arises through cavity-mediated 
cooling 16|, and as we shall see, has an important ef- 
fect on the axial dynamics of atoms in the experiment of 



We now wish to use these potentials and heating rates 
to gain an intuitive understanding of the character of 
atomic motion that we would expect to observe in each 
case. In particular, we are interested in exploring the 
degree to which the atomic motion in the potential can 
be close to conservative motion, or likewise the degree to 
which it could be dominated by diffusion. 

The time scales of relevance to the conservative mo- 
tion may be characterized by the period associated with 
small-amplitude oscillations in the bottom of the axial 
(tq = 1/ fa) and radial (r^ = 1/ ,fr) potential wells. If the 
energy changes only by a small fraction (relative to the 
the total well depth Uq) on this time scale, motion will be 
nearly conservative. Fig. 5 plots the potentials and heat- 
ing rates for the two cases in this new set of scaled units; 
heating rates are expressed as an energy increase per os- 
cillation period, as a fraction of Uq (note as the atom 
heats and explores the anharmonicity of the potential, 
this only lengthens the period of oscillation). Interest- 
ingly, we see a clear qualitative difference in the nature 
of the atomic motional dynamics. For the parameters 
of Hood et al. , in the radial plane spontaneous emission 
only gives small perturbations to the energy over the time 
scale of single orbits, and motion is nearly conservative. 
We note that this low level of diffusion enabled the re- 
constructions of single-atom trajectories in Ref. §, for 
which the small changes in angular momentum could be 
accurately tracked. A quite different regime is found for 
the parameters of Pinkse et al, where the radial atomic 
motion is strongly affected by heating from spontaneous 
emission kicks. Here an average atom gains an energy of 
nearly half the well depth in what would be a radial orbit 
time, adding a large diffusive component to the motion. 
This same scaling shows that the axial heating rate is also 
much more rapid on the scale of the potential in Ref. [|j , 
which suggests that the atom will more quickly escape 
its confinement near an antinode and begin to skip along 
the standing wave. The qualitative understanding of the 
atomic motion gained here is borne out by the simula- 
tions of Refs. j?! and and is explored in more detail 
in the simulations to fo low. 
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B. Simulated transits 



Simulations of the kind described in Sec. II were per- 
formed for the parameters of the two experiments, and 
individual instances of these simulations give insight into 
the dynamics of the motion — for example, the relative 
significance of conservative or dissipative dynamics — and 
the correlation between atomic motion and the cavity 
field state, which is in turn measured by detection of the 
output field. Ensembles of these trajectories provide the 
statistics of the motion described by the Fokker-Planck 
equation [Eq. ^ which may then be used to provide his- 
tograms of transit times to compare to the experimental 
data or to test reconstruction algorithms for the motion. 
In order to approximate the experiment as closely as pos- 
sible, some effort was made to match the detailed exper- 
imental conditions. The two general considerations were 
to reasonably accurately estimate the initial distribution 
of atomic positions and momenta for atoms and to con- 
sider detection noise and bandwidth when simulating the 
feedback switching of the probe laser power. 

For each trajectory in the simulations, initial atomic 
position and momentum values were drawn from a prob- 
ability distribution, which was chosen to correspond to 
the cloud of atoms following laser cooling and then free 
fall or launching by an atomic fountain Q to the cav- 
ity mirrors. In the simulations, all the atoms started in 
a horizontal plane 1 | mode waists above or below the 
center of the cavity mode, where mechanical effects on 
the atom are negligible. Since the MOT from which the 
atoms are falling or rising has dimensions much larger 
than the cavity mode, the initial position in the axial di- 
rection was chosen from a flat distribution over the cavity 
mode, and the initial position along the y axis was also 
chosen from a flat distribution over 1 i mode waists on 
either side of the mode center — this distance could be 
modified but atoms that are far out in the mode radially 
do not typically cause large increases in the cavity trans- 
mission, and therefore do not trigger the feedback. The 
velocity of the atom along the cavity axis is limited by the 
fact that it must not hit one of the mirrors while falling 
toward the cavity, and this was also chosen from a flat dis- 
tribution where the speed was not more than 0.46 cm/s 
for the cavity of Hood et al. 0. Although the two exper- 
iments have rather different geometries, we estimate that 
this consideration leads to a very similar limiting velocity 
for motion along the axis. In the experiment of Pinkse et 
al. 1^, we used 0.4 cm/s. The velocity along the z axis 
was chosen from a Gaussian distribution appropriate to 
the temperature of the MOT (~ 20 ^K) after polariza- 
tion gradient cooling. For Q the velocities in the vertical 
direction were chosen by calculating as appropriate for an 
atom falling freely from the MOT (the MOT is situated 
3.2mm above the mode with a spatial extent of standard 
deviation 0.6mm). Thus atoms arriving at the cavity axis 
have a mean vertical velocity v = 25cm/s. Some of these 
parameters such as the height, size and temperature of 



the initial MOT are not precisely known for the experi- 
ment, so that some consideration of the variation of the 
histograms and other features of the resulting simulations 
has been made although no systematic optimization in 
order to obtain the best agreement has been undertaken. 
In Ref. 1^ the mean initial vertical velocity of atoms en- 
tering the cavity is 20 cm/s. This speed is very much 
less than the mean velocity imparted to the atoms by 
the pushing beam which launches them from the MOT 
25 cm below, and as a result the atoms are all near the 
top of their trajectories. Simple kinematical calculations 
show that the resulting distribution of velocities should 
be rather broad compared to the mean. In the absence of 
more detailed information about the MOT temperature 
and spatial size and the strength of the pushing beam 
we choose the initial vertical velocity distribution to be 
a Gaussian of mean 20 cm/s and standard deviation 10 
cm/s — this leads to a distribution of trapping times with 
a mean that matches the mean reported in Ref. |^ . Each 
trajectory proceeds until the atom is either a greater ra- 
dial distance from the center of the mode than it started 
from, or it has moved sufficiently far in the axial direction 
that it would hit one of the cavity mirrors. 

The detection and triggering are modeled as follows. 
In the parameter range in which the "quasiclassical" 
model is valid, the cavity fleld comes to equilibrium with 
the atomic position on a time scale much faster than 
the atomic motion itself, and thus the light transmitted 
through the cavity (over bandwidths of the order of tens 
to hundreds of kilohertz) is associated with the atomic 
motion. At each point in the simulation the intracavity 
field and intensity expectation values are stored in or- 
der to record for each trajectory a noiseless and infinite- 
bandwidth trace. In practice, experimental traces will 
look like filtered and noisy versions of these traces. As 
an atom enters the cavity mode, a weak driving field is 
present for probing. In order to model the triggering 
step, the field intensity (a^a) or field amplitude modulus 
squared |(a)p is averaged over a time equal to the band- 
width of the detection in the case of heterodyne detection 
as in Ref. 0, or over the time windows in which photo- 
counts are binned in the case of direct photodetection as 
in Ref. . A random number with the appropriate vari- 
ance to represent the shot noise is added and the total 
is compared with some predecided level — if the transmis- 
sion exceeds this level the probe laser beam is increased in 
strength in order to attempt to trap the atom. In the case 
of Ref. [0 the trigger level is |(a)p = 0.32, the averaging 
time is 9 /is, and there is a 2-/is delay between triggering 
and changing the driving laser power. For the experi- 
mental bandwidth of 100 kHz, the appropriate noise has 
standard deviation 0.05 at a transmitted signal of 0.32. 
These parameters are chosen so as to match as closely 
as possible the conditions of the experiment. The same 
procedure is followed for simulations of the parameters 
g. Although the exact triggering protocol is not de- 
scribed there, we assumed that counts over a period of 
10 /iS were used to decide whether or not to trigger and 
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the noise was chosen to be consistent with the reported 
photon count rate of 2 x lO^s^"^ 

Examples of such trajectories are plotted for the pa- 
rameters of Ref. [0 in Fig. 6 and for those of Ref. |8| in 
Fig. 7. The chosen trajectories range in length from the 
experimentally reported mean transit time upward, and 
are chosen because they show typical features of the dy- 
namics in each case. It is clear that the two experiments 
are in quite different parameter regimes, as already indi- 
cated by the relative sizes of the potentials and heating 
rates. 



For the parameters of Ref. [Q, the atoms orbit in a 
radial plane; some have nearly circular and some very 
eccentric orbits. The motion along the axial direction 
is usually well localized near an antinode of the stand- 
ing wave, where the axial heating rate is small. This 
localization occurs because atoms are channeled into the 
antinodes by the weak potential associated with the ini- 
tial probing field, which slowly begins to affect an atom as 
it falls across the mode waist during the detection stage 
of the experiment. However the strong axial heating that 
is present away from the antinodes means that once an 
atom begins to heat axially, it suffers a burst of heat- 
ing (over several hundred microseconds), which leads to 
its loss from the potential well associated with a single 
antinode of the field. Frequently the atom leaves an ax- 
ial potential well when it is radially far from the center 
of the cavity mode, since in this case the axial poten- 
tial becomes weaker. Note that the mean transit time 
in Ref. corresponds to 3.5 radial orbits around the 
center of the cavity mode, so transits with multiple os- 
cillations are frequently observed. In Ref. Q the radial 
oscillation frequency is slower, so an atom of mean tran- 
sit time does not in fact make a complete rotation about 
the mode center. The radial motion in this case is also 
visibly more stochastic in nature, as a result of the rel- 
atively faster spontaneous emission momentum diffusion 
discussed above. 

Another interesting difference between the two param- 
eter regimes is, as suggested in Ref. the relative im- 
portance of atomic motion along the standing wave as 
opposed to oscillations around a single antinode. In the 
case of Ref. [||, long, strongly trapped transits almost 
always involve intervals when an atom is skipping along 
the standing wave, as well as intervals when it is oscillat- 
ing in an individual well. By contrast, for the parameters 
of Ref. 0, only a few percent of trajectories involve skip- 
ping during times in which the atom is trapped, and this 
is usually associated with movement over one or two wells 
with the atom falling back into the adjacent or a nearby 
well. This happens so quickly that it does not affect the 
radial motion in practice, or lead to a detectable signal 
in the output light, so that these rare events of skipping 
do not affect the reconstructions of Ref. ifTl . 
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FIG. 6. Typical trajectories from simulations of the exper- 
iment of Hood et al. as described in the text. The driving 
parameters are Aprobe/2n = —125 MHz and Aac/2n = —47 
MHz, with an empty cavity mean-field strength of m = 0.3 
photons. The trajectories have transit durations of (a) 345 
fj,s, (b) 680 ^s, and (c) 1032 ^s. This is one, two, and three 
mean transit times respectively, (i) The radial trajectory of 
the atom; the z position is plotted against the y position, (ii) 
The y position (dashed line ) and z position (solid line) are 
plotted as a function of time, (iii) The axial position, where 
zero is an antinode of the cavity field, (iv) The noiseless in- 
finite-bandwidth transmission m (solid line) and the radial 
distance from the center of the mode (dashed line). 
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FIG. 7. Typical trajectories from simulations of the exper- 
iment of Pinkse et al. as described in the text. The driving 
parameters are Aprobe/27r = —40 MHz and Aac/iTi ~ —35 
MHz, with empty cavity photon number n = 0.9. The tra- 
jectories have transit durations of (a) 247 ^s, (b) 514 fis, and 
c) 1358 fis. The experimentally reported mean transit time is 
250 /iS. (i) The radial trajectory of the atom; the z position 
is plotted against the y position, (ii) The y position (dashed 
line) and z position (solid line) are plotted as a function of 
time, (iii) The axial position, where zero is the mean axial po- 
sition over the transit, (iv) The noiseless infinite-bandwidth 
transmission n (solid line) and the radial distance from the 
center of the mode (dashed line). 



As noted in Ref. [Q, the axial motion often becomes 
more significant at the end of a transit and as the atom 
is leaving the mode, which leads to atoms skipping a well 
in perhaps as many as one in five cases at the end of the 
transit. We find from the simulations that in Ref. 
the first escape time from an axial potential well for an 
atom initially localized near an antinodc is sufficiently 
short compared to the mean trapping time that skipping 
along the wells almost always takes place. On the other 
hand, the first escape time is of the order of several times 
the mean trapping time for the parameters of Ref. , so 
skipping between standing wells is correspondingly rare. 

It is interesting to note that the friction coefficient for 
the parameters of Pinkse et al. is much more signifi- 
cant than for the experiment of Hood et al, and plays 
an important role in the axial motion of the atom. As 
in the trajectories shown here it is a feature of essen- 
tially every trajectory for the parameters of Ref. Q that 
the atom spends time in potential wells associated with 
several different antinodes of the field. However, we per- 
formed simulations with the sign of the friction coefficient 
reversed, and found that no more than a few percent of 
trajectories were recaptured in a second well after having 
begun to skip along the standing wave. Clearly the dis- 
sipative nature of the motion is an integral feature of the 
dynamics in this regime, and in particular it enables the 
atoms to fall back into axial potential wells after escape 
due to the rapid heating in that dimension. 



IV. SIMULATION RESULTS FOR THE 
EXPERIMENT OF HOOD ET AL. 

Having presented the theoretical basis underlying the 
simulated atom trajectories, in this section we present 
results of these simulations and their comnarison with 
experimental results as reported in Ref. 0]. We gen- 
erate a set of simulated trajectories for the parameters 
{go,l,K) = 27r(110, 2.6, 14.2) MHz with detuning pa- 
rameters Aqc = ujcavity - ujatom = -2?! X 47 MHz and 
^probe = i^probe - ^^atom = "Stt X 125 MHz. In Cor- 
respondence with the experimental protocol, the initial 
pretriggering level of the driving laser gives a 0.05-photon 
mean-field strength in the empty cavity; when this level 
rises to 0.32 photons indicating the presence of an atom, 
we trigger a sixfold increase in the driving strength to 
a trapping level of a 0.3-photon empty-cavity mean field 
strength. A close correspondence between theory and 
experiment is obtained for these results, demonstrating 
the relevance of this theoretical model to the physics of 
the actual experiment. In addition, both theoretical and 
experimental results exhibit features which are relevant 
to building up a picture of the nature of the single-atom, 
single-photon trapping and atomic dynamics, both qual- 
itatively and quantitatively. 

We begin by presenting the qualitative similarity of 
experimental and simulated atom transit signals, as ob- 
served via detection of cavity transmission as a function 
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of time. Fig. 8 shows two sample experimental transits 
[(a) and (b)] and two sample simulated transits [(c) and 
(d)]. For the simulated transits, traces of the correspond- 
ing radial and axial motion are also shown. Transmission 
is shown here as m = | (a) p, as is appropriate for the bal- 
anced heterodyne detection of Rcf. | ^ . In the case of the 
simulated results, the simulated transmission signal has 
been filtered down to the experimental detection band- 
width of 100 kHz, and both technical noise and shot noise 
have been added [Q. The transmission signal thus pro- 
cessed can be seen to lose some of the clarity with which 
it reflects the full atomic dynamics, in comparison to the 
transmission traces of Fig. 6. In particular, the experi- 
mental detection bandwidth is much slower than the time 
scale for axial oscillation in the confining potential, so 
that observed transmission signals are averaged over the 
fast variation in g caused by these axial oscillations. The 
observed maximum transmission should therefore be low- 
ered relative to theoretical predictions, by an amount de- 
pendent on the amplitude of typical axial motion. Thus 
this finite-bandwidth effect allows for an experimental 
estimation of the axial confinement of a typical transit. 
Such a procedure gives an estimate of confinement within 
~ 70 nm of an antinode, in good agreement with simu- 
lation results which suggest typical confinement within 
~ 50 nm. It is important to note that while such tight 
confinement appears typical over the duration of a tra- 
jectory, atoms commonly undergo rapid diffusive heating 
near the end of their confinement lifetime, which leads to 
their escape in a majority of cases. 

FIG. 8. (a) and (b) Examples of atom transits, i.e., cavity 
transmission as a function of time as an atom passes through 
the cavity field for the experiment of Hood et al. Solid traces 
show atoms trapped using the triggering method described, 
with an m ~ 1 photon peak field strength. For comparison, an 
untriggered (untrapped) atom transit is shown in the dashed 
trace. For these traces, the parameters are those of Fig. 3. 
The empty-cavity 0.3-photon mean-field strength is indicated 
by the horizontal dashed line, (c) and (d) Theoretical sim- 
ulation of atom transits for the same Aprobe and Aac- Shot 
noise and technical noise have been added to the transmission 
signals, which have also been filtered to experimental band- 
width. Other traces show the radial (dashed line) and axial 
(solid line) motion of the atom. Motion along x, the stand- 
ing-wave direction, has been multiplied by 10 to be visible on 
the plot. Note that the atom is very tightly confined in x until 
rapid heating in this direction causes the atom to escape. 



A. Trapping lifetimes 

From the entire set of experimental and simulated tra- 
jectories like those of Fig. 8, it is possible to investi- 
gate some quantitative aspects of the trapping dynam- 
ics. First we focus on the trap lifetimes produced by 
the triggered-trapping scheme. Figures 9(a) and 9(b) 
show histograms of experimental transit times for un- 



trapped atoms and for atoms trapped by means of the 
triggered-trapping strategy. Transit durations are deter- 
mined from the experimental data by recording the time 
interval during which the transmission signal is clearly 
distinguishable from the empty-cavity transmission level, 
in the presence of experimental noise. Since the signal-to- 
noise ratio for observing transits depends on the specific 
probe parameters, one must be careful to compare un- 
triggered and triggered transits observed with the same 
detunings and intracavity field strengths. The sole dif- 
ference must be that in the untriggered case, the empty 
cavity field is set at a constant strength so that the atom 
falls through the effective potential, whereas in the trig- 
gered case the field begins at a lower level and is only 
turned up once the atom enters the cavity, thus confin- 
ing the atom. For example, Fig. 8(a) shows sample un- 
triggered (dashed) and triggered (solid) transit signals 
which correspond to one another in this way. 

In Fig. 9 the difference in transit lifetimes between trig- 
gered (b) and untriggered (a) cases is immediately strik- 
ing. For their initial fall velocity oi v = 25 cm/s, atoms 
have a free-fall time of ~ 110 /iS across the cavity waist 
2wo = 2(14.06 /im). As discussed above, the duration of 
observed transits is limited by the signal-to-noise ratio, 
which provides a slightly more restrictive cut on transit 
durations, so the untriggered data set shows a mean du- 
ration of 92 lis. In contrast, when the triggered-trapping 
strategy is employed, the mean trapping lifetime is 340 
fis. The dispersion about the mean likewise changes dras- 
tically from 75 /iS in the untriggered case to 240 /xs in the 
triggered case. These results represent a clear signature 
of the trapping of single atoms with single photons via 
this method. In this setting, atoms have been observed 
to remain trapped in the cavity field for as long as 1.9 
ms. 

The corresponding theoretical histograms are shown 
in Figs. 9(c) and 9(d) for the untriggered and triggered 
cases. The start of the transit is taken to be the time at 
which an atom could be distinguished in the cavity given 
the signal to noise, and the final time is taken to be the 
last point at which the transmission dropped to within 
the noise of the transmission with no atom. This defini- 
tion accounts for the fact that as atoms move out in the 
radial direction the transmission often drops to around 
the free space value, but returns again to some large value 
over the time scales of the atomic motion. These levels 
were chosen to duplicate as closely as possible the proto- 
col for deciding transit times for the experimental data. 

The simulated transit set shows a mean trapping time 
of 96 /iS in the untriggered case and 383 fis in the trig- 
gered case and dispersions of 84 and 240 /iS, respectively. 
This result is in good agreement with the experimental 
results when statistical errors and uncertainties in the ini- 
tial MOT parameters are taken into account. The agree- 
ment between experimental and simulated trap lifetimes, 
in both mean and distribution, gives an indication of the 
validity of the theoretically calculated trapping potential 
and diffusive forces on the atom. The 3.5-fold increase 
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in observed lifetimes due to trapping is made possible 
by the cavity QED interaction, which allows creation of 
a deep trapping potential without correspondingly large 
diffusion as in the free-space case. 
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FIG. 9. Observed atomic transit durations for untriggered 
and triggered cases, with the parameters of Fig. 3. (a) and 
(c) Experimental data show a mean observation time of 92 /is 
in the untriggered case (a) and 340 fis in the triggered case 

(c) , indicating the significant trapping effect. For comparison, 
the free flight time across the cavity waist is 110 fis. (b) and 

(d) The simulated transit set shows a mean of 96 /is for the 
untriggered case (b) and 383 for the triggered case (d), in 
good agreement with experiment. 



B. Oscillations and radial motion 

We now turn to a more detailed investigation of the dy- 
namics of motion experienced by a trapped atom. As we 
have seen, the transmission signal for a single trapped 
atom exhibits large variations over time which may be 
tentatively identified with atomic motion in the radial 
(Gaussian) dimensions of the cavity field. Thus, for ex- 
ample, the highest transmission occurs when the atom 
passes closest to the cavity axis, p = 0. To determine 
the validity of such an identification, we examine the 
periods of observed oscillation in the transmission sig- 



nal. The calculated effective potential is approximately 
Gaussian in the radial dimension, so a one-dimensional 
conservative-motion model predicts periods as a function 
of oscillation amplitude in this anharmonic effective po- 
tential well. Referring to the sample transits of Fig. 8, 
one does indeed note a trend toward large modulations 
with long periods and smaller modulations with shorter 
periods. To quantify this observation, we plot period P 
versus the amplitude A for individual oscillations, where 
A = 2[(iJi + H2)/2 - H,]/{H^ + i?2), with {H,,H2,H,} 
as indicated in Fig. 8. 
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FIG. 10. (a) For experimental transmission data of Hood 
et al., the modulation period is shown as a function of ampli- 
tude. If modulations in transmission are tentatively identified 
with radial atomic motion, their expected period is half that of 
the radial motion. The solid curve gives calculated period vs 
amplitude based on this assumption and on one-dimensional 
motion in the effective potential U{p,0) of Fig. 3. (b) Corre- 
sponding plot for simulated transmission data. Points with 
lowest underlying atomic angular momentum are plotted with 
circles; separation by angular momentum reflects the conser- 
vative nature of atomic dynamics on time scales comparable 
to a radial period. 
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Figure 10(a) shows the experimental data plotted 
along with the calculated curve for one-dimensional mo- 
tion in the effective potential U{p, 0) (see Fig. 3), for the 
same parameters as Fig. 8. (This is a different data set 
from that presented in Fig. 4. of Ref ||^.) Note that since 
an atom approaches the cavity axis p = twice over 
the course of one orbital period, the predicted period 
for oscillations in the transmission signal is half the pe- 
riod of the underlying atomic motion. Experimental data 
clearly map out this calculated curve for radial atomic 
motion, demonstrating that oscillations in the observed 
cavity transmission do indeed reflect radial position of an 
atom as it varies over time within the trap. The agree- 
ment also indicates the quantitative correctness of the 
theoretical model for the radial potential depth and spa- 
tial profile. Note that the comparison is absolute with 
no adjustable parameters. 

The same analysis may be performed for transmission 
oscillations in the set of simulated transits, yielding the 
plot of Fig. 10(b). This plot again shows agreement with 
the calculated curve, with some spread away from the 
line. For simulated transits, it is possible to turn to the 
underlying atomic position record to determine an an- 
gular momentum for the atom during a given oscilla- 
tion. Thus the oscillation data of Fig. 10(b) are plot- 
ted by atomic angular momentum, where lower angular 
momentum data points are shown with circles. A sep- 
aration by angular momentum is clearly evident, with 
lower angular momentum points most closely following 
the calculated one-dimensional (and thus zero angular 
momentum) curve. This separation, while it may seem 
expected, is in fact a non-trivial indication that angular 
momentum is a valid quantity for the atomic motion over 
the course of an oscillation period. Since the atomic mo- 
tion is not in fact conservative, but is also influenced by 
random (diffusive) forces, a separation by angular mo- 
mentum can only be expected to occur if the effect of 
diffusive forces is sufficiently small over the time scale of 
an orbit in the conservative potential. The plots of Fig. 
3 provide an initial indication that this is indeed the case 
for these parameters, and this idea is borne out by the 
current investigation. Confidence in the relatively small 
effect of diffusion over a single orbital period is crucial in 
the reconstruction of two-dimensional atomic trajectories 
as in Ref. 0. 



V. SIMULATION RESULTS FOR THE 
EXPERIMENT OF PINKSE ET AL. 

Having provided a validation of our capabilities for nu- 
merical simulation by way of the results of Sec. IV, we 
next apply this formalism to the experiment reported in 
Ref. [||. At the outset, we note that the various approx- 
imations discussed in Sec. II related to the derivation 
of this quasiclassical model are satisfied to a better de- 
gree for this experiment than for the experiment of Ref. 
m. Hence we expect that the correspondence between 



the simulations and experiment should be at least of the 
quality as in the preceding section. 

Our starting point is the generation of a large set of 
simulated trajectories for the parameters reported in Ref. 
B, namely, (go,7,'«) = 27r(16,3, 1.4) MHz with detuning 



^cavity ^atom 



^27r X 35 MHz and 



parameters A^c 
Apro&e = ujprobe - i^atom = "Stt X 40 MHz. The initial 
pretriggering level of the driving laser gives a 0.15-photon 
mean intensity in the empty cavity; when this level rises 
to 0.85 photons, indicating the presence of an atom, we 
trigger an increase in the driving strength to a trapping 
level of 0.9-photon empty-cavity intensity. These cri- 
teria are intended to follow the parameters indicated in 
Figs. 2 and 3 of Ref. Q. Note that for the cavity geom- 
etry of this experiment, the time for an atom to transit 
freely through the cavity mode in the absence of any light 
forces is tq = 2wo/v — 290 ^s, where as before we take 
twice the cavity waist wq as a measure of the transverse 
dimension of the cavity. 



A. Histograms of transit durations 



From the set of such simulated trajectories (~ 400 in 
this particular case), we can construct histograms for the 
number of events as a function of total transit signal du- 
ration. Following the experimental protocol of Ref. 
which employed photon counting, we base this analysis 
upon the intracavity _photon number n = (a^a) rather 



than|(a)p as in Ref. ||7|, although this distinction is not 
critical to any of the following considerations. 
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FIG. 11. Simulated atomic transit durations for untrig- 
gered and triggered cases, with the parameters of Pinkse et 
al., as in Fig. 7. (a) The untriggered transit set shows a 
mean observation time of 160 fis. (b) The triggered transit 
set shows mean duration 280 /xs, in good agreement with the 
experimentally quoted mean of 250±50ps. For comparison, 
free-fall time across the cavity waist is 290 ps. 
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The resulting histograms for the experiment of Ref . ||] 
are displayed in Fig. 11 for the cases of untriggered and 
triggered trajectories. As in the discussion of Fig. 9, the 
external drive strengths are set to be equal for this com- 
parison to provide equal detectability for an atom passing 
through the cavity mode. Detection with lower exter- 
nal drive strength gives a lower signal-to-noise ratio for 
atom detection, which results in detected transit dura- 
tions much shorter than the actual passage time through 
the cavity (which is of order tq = 2wo jv), as for example 
in Fig. 2(a) of Ref. g. 

In support of the validity of our simulations for the 
experiment of Pinkse et al. (including the initial atomic 
velocity and position distribution and the triggering con- 
ditions), note that the mean of 280 /^s for the histogram 
in the triggered case of Fig. 11(b) corresponds quite well 
with that quoted in Ref. B, namely, foxp — 250 ± 50 /is. 
Further, the histograms in Fig. 11 exhibit an extension of 
the mean transit duration from 160 /is for the case of no 
triggering in (a) to 280 /iS with triggering in (b), in sup- 
port of the claim of trapping in Ref. . The dispersion 
of events around the mean is quite large in both cases, 
161 /U s in the untriggered set and 282 ^s in the triggered 
set. The increase in the mean is largely associated with 
an increase in the number of events in the range 200-300 
/iS, as well as in the number of rare events much longer 
than the mean duration. Once again we note that the 
dissipative nature of the dynamics plays a crucial role in 
the observed motion for the experiment of Pinkse et al. 
A histogram of transit durations calculated with the sign 
of the friction coefficient reversed has a lower mean than 
that of transits with no triggering. 

However, it is certainly worth noting that the observed 
"average trapping time" fexp = 250 ± 50 ^s quoted in 
Ref. H, as well as the corresponding mean time from 
our simulations, are smaller than the time tq = 290 /is 
for an atom to transit freely through the cavity mode. 
Additionally, even in the case of no triggering, there is 
already a significant number of events with similar long 
duration to those in (b) with triggering. Such events 
arise from the relatively large contribution of diffusion- 
driven fluctuations whereby an atom randomly loses a 
large fraction of its initial kinetic energy as it enters the 
cavity. That such fluctuations play a critical role should 
already be clear from the plots of the confining potentials 
and diffusion coefficients in Fig. 4. 



B. Radial motion 

Trapping dynamics can also be explored if atomic oscil- 
lation in the trapping potential can be directly observed. 
Certainly the observations presented in Fig. 10 make this 
case for the experiment of Ref. |Q, with the observed 
oscillation frequencies found to be in good quantitative 
agreement with those computed directly from the anhar- 
monic potential without adjustable parameters and with 



the results of the numerical simulations. 

Towards the goal of constructing a similar plot for the 
parameters of Ref. consider a long-duration transit 
event such as that in Fig. 7(c). Recall that the output 
flux from the cavity is given by the cavity decay rate 2Kd 
into the relevant detection channel times the intracavity 
photon number, or I — 2Kdn — 2Kd{a^a), with then the 
detected count rate found from the overall propagation 
and detection efficiency as R = . Of course, in any 
actual experiment the full information displayed for the 
intracavity photon number n is not available because of 
finite detection efficiencies (^ < 1) and the requirement 
to average over many cavity lifetimes in order to achieve 
an acceptable signal-to-noise ratio (roughly for a time 
such that \/R5t >> 1). 
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FIG. 12. Transmission data for the simulated transit of Fig. 
7(c). The full ideal signal n{t), with infinite bandwidth and 
no degradation due to cavity escape efficiency or subsequent 
system losses, is shown in gray. Slow variations are caused by 
radial motion while fast variations reflect axial motion. The 
black trace results from applying to this ideal signal a low-pass 
filter with cutoff fc ~ 10 kHz intended to optimize the visibil- 
ity of any radial oscillations for frequencies / < 5 kHz, where 

(r) 

/o — 2.6 kHz is the orbital frequency for small- amplitude 
oscillation near the bottom of the radial potential. The re- 
sulting filtered transmission signal shows variations due to 
both radial motion and axial heating. 

Rather than attempt a detailed analysis of such effects 
for the experiment of Ref. § , here we wish to illustrate 
several generic effects that hinder definitive observation 
of radial oscillations in this regime. We therefore take 
the full ideal signal n{t) with no degradation due to cav- 
ity escape efficiency or subsequent system losses (which 
we estimate to be k^/k ^ 0.17 and ^ ~ 0.6for an over- 
all efficiency of 0.11). As shown in Fig. 12, to this ideal 
signal we apply a low-pass filter with cutoff fc = 10 kHz 
intended to optimize the visibility of any radial oscilla- 
tions for frequencies / ^ 5 kHz, where /q*^-* = 2.6 kHz is 
the orbital frequency for small-amplitude oscillation near 
the bottom of the radial potential. As before, recall that 
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a periodic variation in the radial coordinate at frequency 
/ results in a variation in n at 2/. Precisely such a filter- 
ing protocol was implemented for the analysis in Fig. 10, 
there with fc = 25 kHz in correspondence to the larger 
radial oscillation frequencies (/q'"'' = 9.4 kHz for Ref. 0|) 

ii- 

Not surprisingly, the frequent and large bursts of axial 
heating evident for the simulated trajectories of Fig. 7 
result in large variations in the intracavity photon num- 
ber on time scales set by twice the axial oscillation fre- 
quency /o"'* « 430 kHz. While these axial oscillations 
cannot be directly resolved in the detected counting sig- 
nal R{t), their envelope nonetheless leads to variations 
in n{t) and hence R{t) on time scales comparable to that 
associated with radial motion (i.e., 1/2/0*^^), as is appar- 
ent in Fig. 12. Consequently, the low-pass filtering [or, 
equivalently, the time averaging over segments in R{t)] 
that is required experimentally to obtain an acceptable 
signal-to-noise ratio gives rise to observed variations in 
n(t) that can arise from either axial or radial atomic mo- 
tion. In the particular transit shown in Fig. 12, two ap- 
parent variations on time scales ~ 200 /iS are introduced 
by a filtering of the axial motion, whereas the longer 
modulation (~ 600 — /xs duration) does reflect the radial 
position of the atom. This is something of a generic fea- 
ture of the several hundred simulated transits examined; 
shorter-time-scale modulations (< 300 /xs) can reflect ei- 
ther a genuine radial excursion or a filtering of axial mo- 
tion, whereas very long period variations (500-600 ^s) 
are indicative of radial atomic motions. This simply re- 
flects the fact that the bursts of axial motion tend to have 
time scales limited to a few 100 /iS. 

To illustrate these points further, we have constructed 
a plot of period versus normalized amplitude of transmis- 
sion oscillations from our simulations of the experiment 
of Pinkse et al. js), with the result given in Fig. 13. We 
emphasize that the protocol followed is precisely as for 
the analysis that led to Fig. 10(b) for the experiment of 
Hood et al. (see also Fig. 4 of Ref. 0), with the ex- 
ception of the aforementioned reduction in the low-pass 
cutoff frequency. In marked contrast to that case, here 
there is a poor correspondence between the distribution 
of orbital periods from the ensemble of simulated tra- 
jectories and the prediction from the potential obtained 
from Eq. (^. Referring to the discussion of Fig. 12 above, 
we note that about 2/3 of the points in the 100-300 fis 
range result from averaging over axial motion, whereas 
for longer-period (P^300 fis) modulations, 80% of the 
observed points reflect changes in the radial motion, but 
with associated transmission amplitude typically modi- 
fied by the presence of axial motion. The results of Fig. 
13 [which are for the ideal case of h{t) without signal 
degradation due to finite escape and detection efficiency] 
suggest that only in restricted cases can temporal vari- 
ations in R(t) be attributed to radial motion, and not 
instead of (or in addition to) the envelopes of axial heat- 
ing processes. Indeed, such effects are well known in the 



literature, having been previously discussed for the case 
of individual atoms falling through the cavity mode (al- 
beit without triggering or trapping) [ pl| , pT[ . A similar 
conclusion was reached, namely, that axial heating pro- 
cesses contaminate the frequency band associated with 
radial motion, thereby precluding inferences about radial 
motion. For the data presented by Pinkse et al. [^, the 
long (~ 500 /is) time scale of the modulations suggests 
an assignment of these signals to radial motion; however, 
a more detailed characterization of the atom dynamics 
over a larger ensemble of transits should yield this more 
definitively. 

It is also worth noting that the quoted average trap- 
ping time foxp = 250 ± 50 /zs in Ref. j|] is itself less than 

1//q'^^ — 390 /iS, which is shortest time for a full radial 
orbit. Hence any conclusion about motion in the radial 
plane must necessarily be based upon rare events in the 
tail of the histograms of Fig. 11. The rare occurrence of 
these long events is reflected in the small number of data 
points in Fig. 13, which was constructed from the same 
number of simulated transits as Fig. 10(b). 
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FIG. 13. For simulated transmission data corresponding 
to the parameters of Pinkse et al., the modulation period is 
shown as a function of amplitude. If modulations in trans- 
mission are tentatively identified with radial atomic motion, 
their expected period is half that of the radial motion. The 
solid curve gives calculated period vs amplitude based on this 
assumption and on one-dimensional motion in the effective 
potential (7(p, 0) as in Fig. 4. Points with lowest underlying 
atomic angular momentum are plotted with circles. Lack of 
separation by angular momentum reflects the diffusive nature 
of atomic dynamics on time scales comparable to or shorter 
than one radial period. 



C. Axial Motion 

We next turn to analyze motion along the axial di- 
rection, and to the statement of Pinkse et al. |8| that 
Fig. 4 of Ref. [^ "is direct evidence for the atom moving 
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along the cavity axis," as opposed to instances of localiza- 
tion around an antinode for which "hardly any periodic 
structure is visible." In their analysis, Pinkse et al. em- 
ployed a function g'-^-' (e, r, e), whose intention is to pick 
out two-time correlations in intensity, with an enhanced 
signal-to-noise ratio of intensity fluctuations by measur- 
ing coincidences of photon pairs. Here we attempt to in- 
vestigate manifestations of the axial motion independent 
of the details of any specific such function by analyzing 
n{t) directly by way of a windowed fast-Fourier trans- 
form (FFT). More specifically, for each trajectory from 
a large ensemble from our simulations, we apply a FFT 
to the record n(t) with a Banning window centered at 
time U and of total width 25 ^s, with the window then 
offset sequentially to ti+i — ti + 5 fis to cover the whole 
range of a given atomic trajectory. The window width 25 
fis is chosen to be in close correspondence to the record 
length of 20 fis employed by Pinkse et al. Longer win- 
dow widths do not qualitatively change the results of our 
analysis, while a substantially shorter-duration window 
leads to a loss of requisite frequency resolution. 

Two examples from an extended set of such transforms 
are given in Figs. 14 and 15. Parts (a) of each of these 
figures show the mean intracavity photon number n{t), 
the axial coordinate x{t), and a contour plot of the win- 
dowed FFT A/j • (ri) for a single atomic trajectory for the 
parameters of Ref. ||^. Here A/j- (fi) is the windowed FFT 
of n{t) over the entire duration of the trajectory, with 
ti — to + i X 5 iJ,s. Parts (b) of Figs. 14 and 15 compare 
A/t;(ri) for two particular values of ti, namely, at a time 
t flight corresponding to the midst of a flight of the atom 
over several antinodes of the intracavity standing wave 
(i.e., variations in axial coordinate x by several units of 
lA/2) and at a time tiocaiized for which there is apprecia- 
ble heating along the axial direction but for which there 
is no flight (i.e., the atom remains localized within the 
same axial wefl). The times {tfught, tiocaiized) are indi- 
cated by the arrows in the top two panels of parts (a) . 

Perhaps the most striking aspect of the comparison of 
the spectral distributions {Nfiight{^),.N'iocaiized{^)} for 
the cases with and without flight is their remarkable sim- 
ilarity [in (b) of Figs. 14 and 15]. Both display prominent 
peaks near Q,p/2'K = fp ~500-600 kBz, which is in accord 
with the expected frequency for large-amplitude oscilla- 
tion in the axial potential, for which the harmonic fre- 
quency /g°^ ~ 430 kBz (recall that frequency of atomic 
dynamics is half the frequency of the associated varia- 
tions in n(t)). This result is also in accord with that 
from Fig. 4(b) of Pinkse et al., for which their simulation 
leads to 1/rp ~550 kBz for variations in the function g^'^^ 

Bowever, our analysis, as in the comparison of 
{J^fiighti^),Mocaiizedi^)} abovc, indicates that neither 
the observation of a peak in Af{fl) around fip nor of os- 
cillatory structure in ^(''^(e, r, e) around Tp ~ 2'K/Qp is 
sufficient to justify direct evidence for the atom mov- 
ing along the cavity axis. Rather, peaks in A/j . (fi) are 
ubiquitous around frequencies Q.p/2'K ~500-600, and ap- 



pear whether the atom's motion is localized (but heated) 
within a given axial well or whether the atom is in flight 
across several wells. This feature follows from an analysis 
of the full record of n{t) without the deleterious effects of 
finite escape and detection efficiency, or of finite detec- 
tion bandwidth. Such a result suggests that the measure- 
ments of Fig. 4 in Ref ||^ are not in and of themselves 
sufficient to establish unambiguous observation of atomic 
motion across several wells of the cavity field standing 
wave. 
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FIG. 14. (a) Mean intracavity photon number n{t), ax- 
ial position x{t), and a contour plot of the modulus of the 
windowed FFT 7Vt(r2) of n{t) for a simulated transit for the 
parameters of Pinkse et al. (b) At the times indicated in 
(a), \MiocaUzed{^)\ IS plotted corresponding to the arrow at 
tiocaiized = 652 /is (solid curve) and \Mfiight{^)\ correspond- 
ing to the arrow at tfught = 867 /is (dash-dotted curve). 
There are apparently only minor differences between these 
two spectra, which does not support the conclusion about 
axial motion drawn from Fig. 4 in Ref. [8]. 



Our analysis does suggest that it may still be possi- 
ble to distinguish between axial motion confined within 
a well and flight along the cavity axis through a more 
careful quantitative analysis of the respective spectral 
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distributions {Affiight{^),MocaUzed{^)}- With reference 
to Figs. 14 and 15, note that a principal distinction be- 
tween these cases is that in the case of flight there is a 
large decrease of spectral content in the lowest frequency 
components around = 0. This decrease reflects the fact 
that axial skipping causes full-range variation in g, and 
thereby pulls down the time-averaged value of transmis- 
sion n{t). In addition, we note an increase in MfUght{^) 
as compared to NiocaUzed(^) for Fourier components in 
a broad range around D,p/2 and up to fip. The increase 
appears to reflect atomic motion that, during skipping, 
explores the full nonlinear (anharmonic) range of the ax- 
ial potential. These characteristics of the overall spectral 
distributions seem to discriminate more reliably between 
flight and localized heating than does a single-frequency 
peak criterion; they may still offer an avenue for observ- 
ing atomic skips across the standing wave. 
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FIG. 15. (a) Mean intracavity photon number n{t), ax- 
ial position x{t), and a contour plot of the modulus of the 
windowed FFT 7Vt(n) of n{t) for a simulated transit for the 
parameters of Pinkse et al. (b) At the times indicated in 
(a), \MiocaUzed{^)\ is plotted corresponding to the arrow at 
tiocaiized = 673 /is (solid curve) and \Affiight{^)\ correspond- 
ing to the arrow at tfught = 780 ^s (dash-dotted curve). See 
the text for discussion. 



VI. CONCLUSIONS 



A principal objective of this paper has been to inves- 
tigate the extent to which light-induced forces in cavity 
QED are distinct from their free-space counterparts. Our 
perspective has been to seek qualitatively new manifes- 
tations of optical forces at the single-photon level within 
the setting of cavity QED. Note that the importance of a 
quantum character for the relevant fields or phenomena 
is not ensured by the statement that the mean photon 
number n 1, since this is trivially the case in an equiv- 
alent free-space volume for a field of the same intensity 
as that inside the cavity. 

As a starting point, we have presented comparisons 
between the effective potential Uef / (p, x) in cavity QED 
and the corresponding free-space potential, as well as of 
the diffusion coefficients in both contexts (Figs. 3 and 4). 
Perhaps surprisingly, even in a regime of strong coupling 
as in Ref. g, there are only small differences between 
the cavity QED and free-space potentials and diffusion 
coefficients. Note that the comparison of Fig. 4 includes 
"the back action of the atom on the cavity field" 
and yet there are nonetheless no substantive differences 
between the cavity QED and free-space cases for the ex- 
periment of Pinkse ct al. Hence, although the cavity 
QED interactions do bring a substantial advantage for 
atomic detection within the cavity volume, we conclude 
that the claim of trapping an atom with single photons in 
Ref. 1^ involves no new characteristics unique to the cav- 
ity QED environment, with the conservative forces and 
diffusion largely described by the well-known free-space 
theory (Fig. 4). Friction which enhances trapping in this 
regime can be ascribed to cavity-mediated cooling effects 
]l5| , p^ , which are in themselves not uniquely features of 
the quantized-field treatment. However, more analysis is 
required to determine if the observed effects of friction 
do indeed rely on the cavity-field quantization. 

By contrast, for the experiment of Hood et al., a com- 
parison of the free-space theory and its cavity QED coun- 
terpart demonstrates that the usual fiuctuations associ- 
ated with the dipole force along the standing wave are 
suppressed by an order of magnitude. A semiclassical 
treatment of the cavity field yields large diffusions like 
those calculated for the free-space trap. Indeed, if it were 
not for the reduction of heating in the quantized cavity 
QED case, an atom would be trapped for less than the 
period of a single radial orbit before being heated out of 
the well for the parameters of Ref. . Our calculations 
support the conclusion that the suppression in dipole- 
force heating is based upon the Jaynes-Cummings ladder 
of eigenstates for the atom-cavity system, which to our 
knowledge represents qualitatively new physics for opti- 
cal forces at the single-photon level within the setting of 
cavity QED. 

In terms of a more complete analysis, the effective po- 
tential Ueffip, x) and the diffusion coefficient D{p, x) are 
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important ingredients in the quasiclassical theory that 
we have developed for atomic motion in cavity QED. By 
way of detailed, quantitative comparisons with the ex- 
periment of Hood et al. in Sec. IV, we have validated the 
accuracy and utility of our numerical simulations based 
upon the quasiclassical theory. As part of this compar- 
ison, we have demonstrated agreement between experi- 
ment and simulation for histograms of the duration of 
transit events, with mean ft = 340 /iS for the histogram 
in the triggered case of Fig. 9b extended well beyond the 
mean f„ = 92 /^s for the untriggered case. Furthermore, 
ft exceeds the transit time tq = 110 /iS for an atom to 
transit freely through the cavity mode. The simulated 
trajectories of Fig. 6 together with the comparison of Fig. 
10 for the experiment of Hood et al. strongly support the 
conclusion that atomic motion is largely conservative in 
nature, with only smaller contributions from fluctuating 
and velocity-dependent forces. Atomic motion is pre- 
dominantly in radial orbits transverse to the cavity axis. 
The (suppressed) axial heating is important, but only 
towards the end of a given trajectory leading to ejection 
from the trap. Knowledge of the time dependence p{t) 
for the radial coordinate (by way of the detected field 
emerging from the cavity and the solution of the master 
equation) as well of the confining potential U {p, 0) allow 
an algorithm to be implemented for inference of the ac- 
tual atomic trajectory, as demonstrated in Ref. ||^ and 
discussed in greater detail in Ref. ^t\ . 

In the case of Ref. numerical simulations for the 
parameters appropriate to this experiment lead to his- 
tograms with mean 280 ps in the triggered case of Fig. 
11(a) and 160 ps for the untriggered case of Fig. 11(b), 
which should be compared to the time tq — 290 /is for 
an atom to transit freely through the cavity mode in 
this experiment. The simulated transits of Fig. 7 indi- 
cate that atomic motion in this case is dominated by 
diffusion-driven fluctuations in both the radial and ax- 
ial dimensions with friction playing an important role in 
the axial direction. The character of the motion ham- 
pers inference of atomic motion from the record of intra- 
cavity photon number. Axial heating leads to repeated 
large bursts of axial excursions during an atomic transit, 
and hence to large oscillations in the intracavity photon 
number fi{t). The envelopes of these oscillations have 
appreciable Fourier content in the range of interest for 
observation of radial motion, so that there is not an un- 
ambiguous signature for the radial motion in the record 
of n{t) on short time scales, such as those presented in 
Ref. H. Similarly, the result by Pinkse et al. for hopping 
or flights over the antinodes of the cavity standing wave 
is not substantiated by a closer inspection of the Fourier 
content of the relevant signals. As documented in Figs. 
14 and 15, similar signals can be observed for an atom 
localized (but heated) within a single standing-wave well. 
We emphasize that these conclusions concerning the work 
of Ref. 1^ are based upon the analysis of several hundred 
simulated trajectories, apparently well beyond the few 
cases presented in that paper. 



Beyond these comments directed to the prior work of 
Refs.' @,§, we suggest that the capability for numerical 
simulation of the quasiclassical model of atom motion in 
cavity QED should have diverse applications. For ex- 
ample, we are currently applying the simulations to the 
problem of feedback control of atomic motion. Given the 
capability to infer an atomic trajectory in real time, it 
should be possible to apply active feedback to cool the 
motion to the bottom of the effective potential Uef / (r) ■ 
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